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Abstract 

In this article, we explore the kinematics of timelike geodesic congruences in warped five dimen- 
sional bulk spacetimes, with and without thick or thin branes. Beginning with geodesic flows in 
the Randall-Sundrum AdS (Anti de Sitter) geometry without and with branes we find analytical 
expressions for the expansion scalar and comment on the effects of including thin branes on its 
evolution. Later, we move on to congruences in more general warped bulk geometries with a cos- 
mological thick brane and a time-dependent extra dimensional scale. Using analytical expressions 
for the velocity field, we interpret the expansion, shear and rotation (ESR) along the flows, as 
functions of the extra dimensional coordinate. The evolution of a cross-sectional area orthogonal 
to the congruence, as seen from a local observer's point of view, is also shown graphically. Finally, 
the Raychaudhuri and geodesic equations in backgrounds with a thick brane are solved numerically 
in order to figure out the role of initial conditions (prescribed on the ESR) and spacetime curvature 
on the evolution of the ESR. 

PACS numbers: 04.50. +h, 12.10.-g 



* E-mail : suman@cts . iitkgp .ernet . in 
' E-mail: anir@cts . iitkgp . ernet . in 

* E-mail : sayan@cts .iitkgp . ernet . in 



1 



I. INTRODUCTION 



Almost a century ago, in their pioneering research [lj], Kaluza and Klein (KK) proposed 
unification of four dimensional gravity and electromagnetism in a five dimensional gravity 
framework. This proposal raised a fair amount of curiosity, interest and research activity 
among theoretical physicists, on the physics of extra dimensions. The idea of extra spatial 
dimensions appeared in a new incarnation with the advent of Superstring theories 2|, y| . The 
theoretical existence of branes in String theory eventually motivated the hypothesis that we 
may be living on an embedded, timelike submanifold (the brane) of a higher dimensional 
(D > 4) Lorentzian spacetime (warped or unwarped), as assumed in the so-called Arkani- 



-Sundrum (RS) jfj] braneworld models. 



7| on warped braneworlds, published 



Hamed-Dvali-Dimopoulos (ADD) (4j, |5| and Randal 

The seminal work of Randall and Sundrum (RS) |6 
more than a decade ago, refers to the idea of the scale of the extra dimension being spacetime 
dependent, while addressing the issue of stability, in a two-brane scenario. In a single brane 
scenario or from a purely higher dimensional bulk perspective, the space-time dependence of 
the metric function(s) associated with the extra dimensional coordinate(s) basically imply 
that the scale of the extra dimension depends on the on-brane (four dimensional) spacetime 
coordinates. Except for a brief discussion on RS type models, we shall, in this article, mostly 
work with the single brane scenario and a five dimensional bulk. 

In earlier papers js, 9|, geodesies in warped spacetimes have been investigated in detail. 
However, such a study of geodesies alone cannot tell us about the overall local behavior 
of a family of test particles, as observed in the neighbourhood of a freely falling observer. 



This motivates us to study the eyo 
Raychaudhuri equations, in 1955 



ution of geodesic congruences. Since the appearance of 



10| , relativists have discussed and analysed its implications 



in various contexts. In its original incarnation, the Raychaudhuri equations provided the 



basis for the analysis of spacetime singularities in gravitation and cosmology [11] . For 
example, the equation for the expansion and the resulting theorem on geodesic focusing, is 
a crucial ingredient in the proofs of Penrose-Hawking singularity theorems 12|, Il3 |. 

The kinematics of geodesic congruences is characterised by three kinematical quantities: 
isotropic expansion, shear and rotation (henceforth referred as ESR) 14Nl9|. which evolve 
along the flow according to the Raychaudhuri equations. Though mostly quoted and used 
in the context of gravity, these equations by virtue of their geometric nature, have a much 
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wider scope in studying geodesic as well as non-geodesic flows in nature, which may possibly 



arise in diverse contexts (see 



19| for some open issues). Two of the authors here have 



recently used these equations to investigate the kinematics of geodesic flows in stringy black 



hole spacetimes 



20[ and flows on flat and curved deformable media (including elastic and 



viscoelastic media) in detail 2ll. \22\. 



In this article, we attempt to understand the kinematics of geodesic flows in five dimen- 
sional warped bulk spacetimes with and without branes. In Section II we quickly recall the 
background spacetime geometries and geodesies. Section IIIA analyses flows in the RS ge- 
ometry with and without branes. The ESR, as obtained from definitions, for a background 
with a thick brane are discussed in Section IIIB. Numerical solutions of the geodesic and 
Raychaudhuri equations are presented in IIIC. Finally, Section IV contains our conclusions 
and comments. 



II. THE BULK SPACETIMES AND GEODESICS 



The bulk spacetimes we work with are given by the line element [8| (77 is the conformal 
time), 

ds 2 = e^V^-dr/ 2 + dX 2 } + b 2 (r])da 2 . (2.1) 

Table HI shows the chosen functional forms of f(a) (the warp factor), the cosmological [a(r/)] 
and extra dimensional [b(r])] scale factors (following Q]). We mostly prefer to work with 



/(*) 


{ a(v), Kv) } 


— In (cosh ka) 
Decaying warp factor 


Set (A) {277,1 + 1} 
FRW (Radiation dominated) brane 


In (cosh ka) 
Growing warp factor 


Set(B) {^,1-2} 
de Sitter brane 



TABLE I: The four possible combinations of /(cr), a(rj) and b(rj). 



thick branes in order to avoid the discontinuities and delta functions which appear in the 
connection and curvature, for thin branes. However, as we will see later, in some special 
cases (e.g. Einstein spaces) one can indeed solve the Raychaudhuri equations consistently, 
with zero rotation and shear, in the presence of thin branes. It may be noted that these 



models are assumed to represent the evolution of the universe beginning at a finite time 
when both the scales of visible and extra dimension were same. 

The geodesic equations in the spacetime (12.1 1) are almost impossible to solve analytically. 
However, they may be recast as a dynamical system of coupled, ordinary, first order differ- 
ential equations as given by Eqs 3.8-3.11 in [8|. Thus, for simple cases, knowing the first 
integrals one may directly determine the ESR. 

III. RAYCHAUDHURI EQUATION AND ESR VARIABLES 
A. Randall Sundrum warp factor with and without branes 

In Einstein spaces, Rab ~ 9ab an d the equation for the expansion simplifies considerably 



if we assume the shear and rotation as zero. For the RSI scenario [6|, [7H in the absence of 

n 



any brane, the Raychaudhuri equation Q, [15| gives us 

dO 2 A 

where, A is the bulk cosmological constant and M is the five dimensional Planck mass. With 
= Ay (notion of focusing is related to F = 0, F < at finite A), leads to 



simpl 



F + eF = 0, with *=\/^3, (3.2) 

Following 0], A is negative and Eq J3.2l has simple oscillatory solutions such as c\ sin(A;A-|-C2), 
which imply = 4k cot(/cA+C2). Therefore, the nature of focusing or defocusing of geodesies 
in the bulk depends on the initial condition or the value of c^- However, this is the behavior 
of geodesic congruences in the bulk with no branes. If we introduce two 3-branes, the hidden 
brane (with positive tension 24M 3 k) at a = and the visible brane (with equal negative 
tension) at a = tt, Eq J3.2l becomes 

F+[k 2 + 2k{5{a)-5(a-ir)}} F = 0. (3.3) 

Using the following property of Dirac delta function 

^(A)) = E^F7W 1 ' M 
and the first integral of the a geodesic equation 

a = ^C 2 e~ 2 f - 1 with f(a) = -k\a\, (3.5) 



we arrive at 

F + [k 2 + hS(X - A ) - k 2 S(X - K)] F = 0, (3.6) 

where, h = -^=, k 2 = , 2k A = and X n = — VpZ^E i. The solution of 

the above equation is given by [23] 



F(A) = ci sin(A;A + c 2 ) + c 3 e- a|A ~ Ao1 + c^" 1 ^ 1 (3.7) 

where c±, c 2 , c 3 , C4 are arbitrary constants. We then integrate the second order equation for 
F, around the neighborhood of Ao and A T to obtain two algebraic equations for c\ and c 2 . The 
determinant condition for nontrivial solutions of c%, c 2 yields the following transcendental 
equation, 

(2a - £4) (2a + k 2 ) + ^A; 2 e^ 2a(A ^ A()) = 0. (3.8) 
from which a can be obtained numerically. In our case, a ~ k\/2 is a good approximation 



(we have checked this with numerical solutions as well). Fig. 1(a) shows (shaded regions) the 
domain of the parameters c^/ci and c^jc\ (c 2 is taken to be zero), in cases where F = at 
some finite value of A between the two branes. Points on the boundaries of the shaded regions 
correspond to the occurrence of F = at the location of the branes (i.e. A = Ao, A = A,,- 
and the corresponding a through the function <r(A)). In the lightly shaded region, F > 0, 
which implies complete defocusing of geodesies (8 — > +00) at finite A, whereas in the darker 
region, F < 0, i.e. geodesic focusing (6 — > —00) is possible. Note that different values for 
c 2 will result in different parameter space diagrams involving the quantities c^jc\ and c^jc\. 
Eq l3. 71 leads to the modified expansion scalar, due to the presence of the branes, as given by 

= 1 cifc cos(fcA + c 2 ) - c 3 a sgn(A - A )e~ a l A ~ A °l - c 4 a sgn(A - A 7r )e" a l A ' A -l 
ci sin(A;A + c 2 ) + c 3 e- a l A - A °l + c 4 e~ a l A ~ A -l 

Due to the new integration constants, the behavior of geodesic congruences have become 



much richer. One such scenario is shown in Fig. 1(b) where for specific values of these 



parameters, c\ — 1, C3 = 1 and c 4 = —2 which belongs to the dark shaded domain in 



Fig. 1(a), gives rise to occurrence of congruence singularity in between the brane locations. 



B. Calculating ESR from velocity field 

One may derive analytic expressions for the ESR variables directly from the following 
definitions 



151 ] for the expansion O, the shear and the rotation Qab, using the 
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FIG. 1: (a) Phase diagram in the c-^j c\-c^j c\ plane. Points lying in the deeply (lightly) shaded 
region correspond to focusing (defocusing) of geodesies, (b) Evolution of expansion scalar in 
presence (continuous line) and in absence (dotted line) of two 3-branes with c\ = 1, C3 = 1 and 
C4 = —2 which lies in the dark shaded region in (a). Specific values of the parameters are: k = 12, 
C = 2, Xq = 0.087 (1st brane location), = 0.131 (2nd brane location), C2 = and a = 6.9282 
(from Eqi37 



geodesic vector field components obtained in 



e = v A u A , 

^AB = ^ {^BU A - ^aU B ) ■ 



n 



-h AB Q, 



(3.10) 
(3.11) 

(3.12) 



Here, n is the dimension of spacetime and h A s — 9ab ± u a ub is the projection tensor (the 
plus sign is for timelike curves whereas the minus one is for spacelike ones) and u A u A = =Fl. 
Let us now obtain the ESR for some specific cases (for thick branes), where we take one or 
two of the metric functions as constants. 



1. Case 1: a{rj) = b{rf) = constant 

Here, with only a non-constant /(cr) present in 12.11 the velocity field components for 
timelike geodesies are 

u a = C a e~ 2f where a = 0,1, 2, 3 (3.13) 

3 

M 4 = VC 2 e- 2 / - 1 where C 2 = C 2 - C i , ( 3 - 14 ) 

i=l 
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where the C Q 's are integration constants, which are constrained by the fact that u A has to 
be real valued. 

According to the definitions Eq J3.1(JI - Eq J3.12} the expansion scalar and the other ESR 
variables for a congruence of timelike geodesies, written as functions of a, give us, 

= 4 f, (3.15) 

VC 2 e- 2 f -1 

Ofji -4/ £12 

E 2 = S AB S AB = 4(c2e _ 2/ _ 1} and n AB = 0VA,B. (3.16) 

For the chosen growing and decaying warp factors respectively, the expansion scalar becomes 

3C 2 sech 2 cr - 4 , , ^ 3 C 2 cosh 2 cr - 4 
0+ = — tanh a and 0_ = tanner. (3-17) 

VC 2 sech 2 a - 1 VC 2 cosh 2 a - 1 

It can be easily seen that, for & > and C > 1, 0± — >■ — oo as cr increases. Therefore a 

geodesic congruence singularity arises in both the cases. With 9+, if C > v/4/3, initially 

the expansion remains positive but eventually geodesies meet exactly at the boundary of 

the accessible domain along the extra dimension because the velocity component along cr, 

m 4 , which appears in the denominator in the expression for G, vanishes at that point (and 

also changes sign). With 0_, the geodesic congruence singularity appears as a —> oo. For 

cV < an exactly similar behavior is obtained. From the nature of corresponding geodesies 



(see |8|), we note that + experiences a finite time singularity (i.e. finite A as well as finite 
cr) but 0_ becomes singular at finite A but rj, a — > oo. 

To understand how geodesic congruences behave in all the abovementioned scenarios, let 
us consider the evolution of the projections of the cross-sectional area orthogonal to the flow 
lines, of a congruence of four geodesies, on different two dimensional surfaces. This is done 
by numerically solving 

u c V c Bab = -B AC B C B - R AC bdu c u d . (3.18) 
i%u B = B A B e (3.19) 

along with the geodesic equations. Here, Bab = ^bUa is the gradient of velocity field and 
£ A represents the separation between two neighboring geodesies. Eq J3. 191 is essentially the 
evolution equation for the deviation vector. To see the evolution from a local observer's 
viewpoint, we have to express the tensorial quantities in the frame basis. The metric tensor 
in coordinate basis and frame basis are related as 

9 AB = e A a e B bV ab , (3.20) 



where the vierbein field, e£, has two indices, 11 A" labels the general spacetime coordinate 
(w.r.t. the coordinate basis) and "a" labels the local Lorentz spacetime or local laboratory 
coordinates (w.r.t. the frame basis). The tensorial components in these two bases are related 
as 

t A = e A a t a . (3.21) 

In the frame basis, we set the initial conditions such that, at A = 0, the projected area has 
the shape of a square in the -^/gu^-y/g^C, 2 plane or in the y/gn$, 1 --Jgu^ 4 plane where £'s are 
essentially solutions of Eq l3.19l Fig J2] shows the evolution of the area elements as A increases. 
The origin represents the location of the observer. We have chosen initial conditions such 
that all components of the rotation vanish. 
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(a) Evolution of the projected square element on y / <?i7£; 1 -y / 522£ 2 plane 
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(b) Evolution of the projected square element on y/gui, 1 - v/sii^C 4 plane 
FIG. 2: Evolution of different 2D surface elements in Case 1 in presence of a growing warp factor 

Figj2] shows how a projected 2D square element evolves, in a static bulk, in presence 



8 



of a growing warp factor. In Fig. 2(a) the initial square area (at A = 0) in the y^IT^ 1 - 
Jlfn^ 2 plane expands and distorts slightly and converges on a parallelogram at A ~ 3.14. 



In Fig 2(b) , however the area in the -s/gii^-y/glAC, 4 P iane shrinks, distorts and converges on 
the y/gn!; 1 axis at A ~ 3.14, clearly suggesting focusing along the extra dimension. The 
effect of shear is evident from the evolution of the shape of the area. 
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(a) Evolution of the projected square element on y/giiS, -J()2i£, plane 
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(b) Evolution of the projected square element on y/gn^-y/g^ 4 plane 
FIG. 3: Evolution of different 2D surface elements in Case 1 in presence of a decaying warp factor. 



On the other hand, Figj3] corresponds to the scenario where the warp factor is of decaying 
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type. In Fig. 3(a) shrinking of the area element is quite prominent whereas in Fig, 3(b) 



it is 



not so. However, in both the figures, the square area becomes more and more parallelogram 
shaped and eventually converge on a line y/gu^ 1 oc ^/g-ri^ 2 or y^Ii^ 1 oc y^iJ^ 4 as A — > oo. 

The nature of evolution of the square elements is therefore a distinguishing feature be- 
tween bulk universes with growing and decaying warp factors. It is worth mentioning here 



that, to a brane based observer, only the evolution depicted in Fig. 2(a) or Fig. 3(a) will be 



visible whereas congruence singularities realised in Fig. 2(b) or Fig, 3(b) will remain unno- 
ticed. In a way, therefore, one can find the nature of warping as well as the existence of 
extra dimensions from the evolution pattern of cross-sectional area elements. 



2. Case 2: f{o)= constant 



In this case, with only non-constant a(r]), b(r]) present in 12.11 the velocity field compo- 
nents for timelike geodesies are 



u 



u 



a?' 



u 



9± 
b 2 ' 



\ 



1 



i=l 



a 



a 2 b 2 



(3.22) 



We can calculate the expansion scalar which turns out to be, 



6 = 

u 



Cfa 



3C 2 a 
ab 2 



b v-Cffo 
b + ^ a 2 b 



(3.23) 



As done before one can also find the components of the shear tensor (not shown here). The 
rotation tensor components are zero, as is evident from the velocity field. 

In case of Set (A) [as given in Table H], the geodesies become parallel as 77 — >■ oo. This 
is because, with a FRW (radiation dominated) brane, the expansion of the universe itself 
slows down with increasing 77 (it is worth mentioning here that as the cosmological evo- 
lution is assumed to begin at a finite 77, the past singularity will not appear here since it 
corresponds to 77 = and falls outside the domain of 77 for the models considered in this 
article). For Set(B) [see Table H], the geodesies spread apart at an ever increasing rate with 
increasing 77 - this is due to very rapid (exponential in real time) expansion of the universe. 
In both these cases geodesic focusing is not achieved. On the other hand, as noted in the 
previous subsection, geodesic singularities are unavoidable when a non-constant warp factor 
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is assumed. Therefore, it should be interesting to see how these two apparently opposite 
features compete with each other when we consider the full braneworld geometry with all 
three non-constant metric functions present. 

We end this subsection by trying to figure out the individual effect of the dynamic nature 
of b(rj). Assuming a(rj) as constant in Eq J3.23l we get, 

v / (l + E,C?)6 2 + Q 

In our models b is negative and as t] increases it tends to zero. For the b(r]) of Set (A), O 
tends to zero i.e. geodesies become parallel while for Set(B) it converges to a finite value 
which implies that the geodesies keep moving away from each other at an approximately 
steady rate. So, with only a b(rj), in both the above Sets (A) and (B), geodesic focusing 
never happens. It is interesting to note that even if b(r]) — > i.e. size of the extra dimension 
becomes singular the expansion scalar remains finite as long as b is finite. Therefore we 
expect b(r]) to play a role only in introducing a scaling effect. This will become clearer in 
the next section, when we consider the general scenario where all the three non-constant 
metric functions are considered. 

One may ask-what about analytic expressions for the kinematic variables in the general 
case? Analytic expressions for the first integrals of the geodesic equations with non-constant 
f(a) and a{rf) can be found easily. However, the ESR as obtained from the definitions 
(Eq J3.10tl3.12l) are functions of a, r] and cannot be reduced to explicit functions of A alone. 
This happens because we do not know how <r(A) and 77(A) are related analytically. Thus, 
the above two subcases are the only ones where one can find useful, closed-form analytic 
expressions for the ESR, directly from the velocity field. 

C. Numerical solutions 

Here, we numerically solve Eq J3 . 1 8 1 simultaneously with the geodesic equation for different 
combinations of the metric functions, in order to understand the interplay amongst all 

nn 

the terms appearing in the Raychaudhuri equations pjj, [15|] . The effect of the shear and 
rotation on the expansion are mutually opposite. Thus, the rotation may play a role in 
avoiding/delaying congruence singularities. On the other hand, the curvature term Rabu a u b 
also has a significant effect on ESR profiles through its large positive or negative value at a 

11 



given spacetime point. 

We have analysed each case for two types of initial conditions - one with zero initial ro- 
tation and one with very high initial rotation, keeping initial and Sab as zero. Non-zero 
initial values for Qab are chosen such that £Iabu b = = u a Qab at A = 0. Initial velocities 
for cases involving Set(A) scale factors are taken as {0.728,0.1,0.1,0.1,0.5}, whereas for 
Set(B), it is assumed as {1.1314,0.1,0.1,0.1,0.5} so that the timelike constraint is satis- 
fied. In FigJH we present the evolution of the expansion scalar for the different scenarios 
mentioned in Table HI 
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(c) With /(cr) = log(cosh cr), 
a{rf) = 1/(1 - rf) and 6(77) = 1-77/2 
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FIG. 4: Nature of expansion scalar for different cosmological braneworld scenarios with two different 
sets of initial conditions. The subscripts (h) correspond to zero (high) initial rotation. 



In the presence of a growing warp factor, a FRW (radiation dominated) brane and an 
asymptotically static extra dimension, geodesic congruences without any initial rotation, 



expand slowly at first but later become focused at a finite A, as observed in Fig. 4(a) The 
shear (S 2 ) is found to grow indefinitely (figure not shown). When high initial rotation is 
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introduced, the congruence diverges for a very short while but eventually the geodesies get 
focused again at another finite but larger value of A. Though rotation increases at late times, 
it is always dominated by shear, which grows even faster. Initial rotation only succeeds in 
delaying the focusing. The curvature term, as we note, is not an important factor here. 



In contrast, Fig. 4(b) shows that, in the presence of a decaying warp factor, a FRW 
(radiation dominated) brane and an asymptotically static extra dimension, the geodesies, 
without any initial rotation, do come closer to each other monotonically, but focus only at 
a — > oo. As before, the shear grows indefinitely. When high initial rotation is introduced, 
the congruence expands initially but eventually geodesies tend to focus again asymptotically. 
The curvature term becomes large positive valued and thus always seems to help in the 
occurrence of congruence singularities. 



Fig, 4(c) represents almost the same behavior as seen in Fig. 4 (a) However, here the large 
initial rotation decays down very quickly. The physical reason behind this is the very fast 
expansion of spacetime that smears out all the initial rotation components. 

On the other hand, in the presence of a decaying warp factor, a de Sitter brane with an 



asymptotically static extra dimension, in Fig. 4(d) , we have an example where geodesies are 
defocused irrespective of the initial rotation. Even though it seems that at late times shear 
totally dominates over rotation, it is the curvature term in the Raychaudhuri equation that 
becomes dominant and causes the defocusing. Remember that, with a(rj) = b(rj) = constant, 
congruence singularity was inevitable. On the other hand, we have observed that initially the 
curvature term is very small. This implies that, in this case, a high enough initial negative 
expansion should lead to a congruence singularity at a finite A (before the curvature term 
becomes dominant). This behavior has been checked with an initial expansion, O = —30 
(figure not shown). 

From the results in FigJlJ one can draw some general conclusions about the nature of the 
ESR variables. Congruence singularity is inevitable (irrespective of the initial rotation) in 



the cases addressed in Fig 4(a) and Fig,4(c) This is because, in the presence of growing 
warp factor, the geodesies have a turning point in the extra dimension. This forces the 



congruence singularity to occur. Fig, 4(b) represents a case where the geodesies are not 
bounded. Even in this case, a high initial expansion cannot make the congruence to diverge. 
On the other hand, the geodesies are divergent when f(a) = — log(cosh a), a(rj) = 1/(1 — r/) 
and b{rj) = 1 — rj/2 [Fig. 4(d)] , which corresponds to a negatively warped and exponentially 
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expanding brane. Here, it is interesting to note that the defocusing occurs because of the 
curvature term which becomes dominant and large negative, and has nothing to do with 
the initial rotation. Therefore, one can say that this defocusing is purely an effect of the 
spacetime geometry 
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(a) Evolution of 2D surface elements for /(tr) = log(cosh cr), a(r]) = 2i] 
and b{rf) = 1 + 1/r] with 6 = E 2 = and O 2 - 54 at A = 0. 
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(b) Evolution of 2D surface elements for /(cr) = — log(cosha), 
a(r?) = 1/(1 - rf) and 6(77) = 1 - r//2 with 6 = S 2 = and ft 2 - 1080 at A = 0. 

FIG. 5: Evolution of the projected square elements on Jgu^-Jga^ plane. 



Let us now look at the evolution of a congruence of geodesies from the local observer's 
point of view and see how the ESR profiles, plotted in Fig. HJ are realised. As done earlier, 
here also we have plotted the evolution of square elements, orthogonal to the congruence of 
the timelike geodesies, projected on y/gu^-y/gii^, 4 planes. These plots provide a different 

14 



perspective of the evolution since they involve different shear and rotation tensor compo- 
nents. We have labeled one point of each area element (the one in the second quadrant, 
initially) as "a" . Following the location of this labeled point on future quadrilaterals (with 
increasing A) one can see the effect of rotation. 



Fig, 5 (a) corresponds to the evolution of the dashed curve in Fig. 4(a) (i.e. in presence of 



a growing warp factor). In Fig 5(a) , shrinking of the area element after an initial expansion, 
increase in the amount of shear after a decrease in the middle and a quick decrease in rotation 
are more clearly visible. 



Fig, 5(b) represents the evolution of the dashed curve in Fig, 4(d) (i.e. in presence of a 



decaying warp factor). The area of the quadrilateral keeps on increasing with A. Effect 



of high initial rotation is prominent, so is its rapid decrease. Effect of shear in Fig 5(b) 



matches with the profile of Fig, 4(d) As mentioned earlier, these qualitatively different 



ESR evolutions are in one to one correspondence with different models. Therefore, loosely 
speaking, the pictorial visualisation provides pointers toward possible verification of those 
models, though much more needs to be done in order to arrive at explicit verifiable signatures. 



Warp factor 

p2/(<7) 



Constant a(r]), b(rj) 
(Analytical results) 



a(rj) = 2rj, b(rj) = 1 + I/77 
(FRW Radiation dominated) 



a(rj) = 1/(1 - n), b(r]) = 1 - V /2 
(de Sitter) 



Growing 



Congruence singularity 
at finite a 



Congruence singularity 
at finite a [Fig 4(a)| 



Congruence singularity 
at finite a [Figjgg] 



Decaying 



Congruence singularity 

at a — > 00 



Congruence singularity 
at cr ->■ 00 [Fig ]4(b)| 



Defocusing at a — >■ 00 [Fig 4(d)] 
or congruence singularity at 
finite a for large, -ve initial G 



Constant 



No congruence singularity 



No congruence singularity 



TABLE II: Summary of behaviour of geodesic congruences for different metric coefficients 



IV. DISCUSSION 

We began by considering geodesic flows in the RS background. Without branes congru- 
ence singularities always occur whereas with two branes, the expansion profile indicates how 
focusing/defocusing can occur in the spacetime between the branes. 
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Later, we analyse geodesic flows in a bulk geometry with a thick brane. Using first 
integrals of geodesic motion, analytic expressions for the kinematic variables are obtained. 
We show how differences arise as we change the warp factor from growing to decaying or 
when we do not have any warping but retain the time-evolving cosmological and extra 
dimensional scales. 

Further, we numerically solve the Raychaudhuri and geodesic equations to obtain the 
expansion, shear, rotation and demonstrate the role of initial conditions on their evolution. 
With a(rj) = rj, a growing warp factor leads to a finite 77 (and a) congruence singularity 
whereas for a decaying warp factor, geodesies are focused at 77 (and a) — > 00. For a de Sitter 
universe, a decaying warp factor may fail to focus the geodesies (due to the large negativity 
of the curvature term in the Raychaudhuri equation), though this is not the case with a 
growing warp factor. When the curvature effect is relatively small, a congruence singularity 
can still arise but for high enough negative initial expansion. The effect of initial rotation, 
on the ESR profiles, especially the expansion scalar, is found to be quantitative. Focusing 
without and with initial rotation yields similar results, though with large initial rotation, 
geodesies tend to spread for a while (focusing at larger A value). All our conclusions are 
summarised in Table [TTJ 

For a visual perspective, we have shown snap-shots of the evolution of a square element 
orthogonal to a geodesic congruence, from a local observer's point of view. The evolution 
of the expansion, shear and rotation, along the congruences become more explicit through 
these figures. 

The effect of b{rf) seems to be largely quantitative since, in our models, as A evolves b{rj) 
tends to a static value with a deceleration. 

It may be asked-what relevance, if any, does a congruence singularity have in the context 
of realistic scenarios? After all, congruence singularities are not real spacetime singularities 
where curvatures diverge. Here, we are tempted to draw an analogy from null geodesic con- 
gruences, for which congruence singularities are nothing but the well-studied caustics where 
optical intensities get magnified immensely. Similarly, in the case of timelike geodesies, we 
may end up with accretion-like effects resulting out of matter accumulation in the neighbor- 
hood of a point. For instance, our visualisation analyses do show how the square elements 
change shape, get rotated because of variations in the metric functions. We may contemplate 
such accretion effects for flows around brane-world black holes 24| and in such situations, 



1(3 



it will become necessary to pursue a line of thought very similar to what we have followed 
in this article. 

Acknowledgments 

SG thanks Council for Scientific & Industrial Research (CSIR), India for providing finan- 
cial support and Centre for Theoretical Studies, IIT Kharagpur, India for allowing him to 
use its research facilities. 



[1] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Berlin. (Math. Phys.), 966-972 (1921); O. Klein, 
Z. Phys. 37 (1926) 895 

[2] M. S. Green, J. H. Schwarz and E. Witten, Superstring theory (Cambridge University Press, 

UK, 1987); J. Polchinski, String theory (Cambridge University Press, UK, 1997). 
[3] I. Antoniadis, Phys. Lett. B 246 (1990) 377. 

[4] N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 429, 263-272 (1998); N. Arkani- 

Hamed, S. Dimopoulos, G. Dvali, Phys. Rev. D 59, 086004 (1999). 
[5] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436 (1998) 

257 [arXiv:hep-ph/9804398] . 
[6] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999). 
[7] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 
[8] S. Ghosh, S. Kar and H. Nandan, Phys. Rev. D 82, 024040 (2010). 
[9] S. S. Seahra, Phys. Rev. D68, 104027 (2003). 
[10] A. K. Raychaudhuri, Phys. Rev. 98, 1123 (1955); 

[11] S. W. Hawking and G. F. R. Ellis, The large scale structure of spacetime (Cambridge University 

Press, Cambridge, UK, 1973). 
[12] R. Penrose, Phys. Rev. Lett. 14, 57 (1965). 

[13] S. W. Hawking, Phys. Rev. Lett. 15, 689 (1965); ibid. 17, 444 (1966). 

[14] E. Poisson, A relativists' toolkit: the mathematics of black hole mechanics (Cambridge Uni- 
versity Press, UK, 2004). 
[15] R. M. Wald, General Relativity (University of Chicago Press, Chicago, USA, 1984). 

17 



[16] P. S. Joshi, Global aspects in gravitation and cosmology (Oxford University Press, Oxford, UK, 
1997). 

[17] G. F. R. Ellis in General Relativity and Cosmology, International School of Physics, Enrico 

Fermi-Course XL VII (Academic Press, New York, 1971). 
[18] I. Ciufolini and J. A. Wheeler, Gravitation and inertia (Princeton University Press, Princeton, 

USA, 1995). 

[19] S. Kar and S. SenGupta, Pramana 69, 49 (2007) and references therein; S. Kar, Resonance, 
Journal of Science Education 13, 319 (2008). 

[20] A. Dasgupta, H. Nandan and S. Kar, Phys. Rev. D 79 124004 (2009). 

[21] A. Dasgupta, H. Nandan and S. Kar, Annals of Physics 323, 1621 (2008). 

[22] A. Dasgupta, H. Nandan and S. Kar, Int. J. of Geom. Meth. Mod. Phys. 6(4) (2009). 

[23] T. C. Scott and R. B. Mann, General Relativity and Quantum Mechanics: Towards a Gener- 
alization of the Lambert W Function, AAECC (Applicable Algebra in Engineering, Commu- 
nication and Computing), vol. 16, no. 6, (2006) |arXiv:math-ph/060701l] y2] . 

[24] C. S. J. Pun, Z. Kovacs and T. Harko, Phys. Rev. D 78 (2008) 084015 



18 



